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Abstract 



Two-spin correlations generated by interactions which decay with distance 
r as r~'^~" with — 1 < ex < are calculated for periodic Ising chains of 
length L. Mean-field theory indicates that the correlations, C{r,L), diminish 
in the thermodynamic limit L ^ oo, but they contain a singular structure 
for r/L ^ which can be observed by introducing magnified correlations, 
LC{r, L)—J2r C{r, L). The magnified correlations are shown to have a scaling 
form <I>(r/L) and the singular structure of ^(x) for x ^ is found to be 
the same at all temperatures including the critical point. These conclusions 
are supported by the results of Monte Carlo simulations for systems with 
a = —0.50 and —0.25 both at the critical temperature T = Tc and at T = 2Tc. 
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I. INTRODUCTION 



The Ising model in d = 1 dimension with long-range interactions decaying with distance 
r as J(r) ~ f^^^'^ has been much studied. Dyson |l| showed that there is no phase transition 
for cr > 1, while there is long-range order at low temperatures when the interactions are 
ferromagnetic and cr < 1. The borderline case of cr = 1 was studied by several authors 0, 
who found a one-dimensional version of the Kosterlitz-Thouless transition ||^ in the system. 
A variety of methods have been used to show that the transition is of second-order 
with continuously varying critical exponents in the region ^ < cr < 1. For < cr < |, it has 
been found that the critical exponents are mean field like for thermodynamic quantities, but 
the spin-spin correlations at the critical temperature 

oc — 4- = -\ (1) 

have a non-classical critical exponent = 2 — a . The case cr = | is somewhat special in 
that the mean field exponents are modified by logarithmic corrections in analogy to the 
situation with the short range Ising model for d = 4, or the d = 3 Ising model with dipolar 
interactions [|1^ . 



The case cr < does not appear to have been investigated except for cr = — 1. The 
reason, perhaps, is the fact that the interaction potential, J{r), becomes nonintegrable for 
cr < and then the energy of the system is non-extensive. The textbook example of cr = —1 
has attracted much attention since it corresponds to an Ising model in which every spin 
interacts equally with every other spin. The problem of non-extensive energy is solved by 
rescaling the coupling strength by the number of spins and the result is a mean-field exact 
solution. Spatial correlations have no meaning in this system. 

For cr < and cr 7^ — 1 one can still consider spatial correlations and this is what we 
shall do in this paper. The reason for doing it is partly the fact that there are interactions 
in nature which are non-integrable. The gravitational interaction is a prominent example 
but there are many systems which contain topological defects with non-integrable effective 
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interactions between the defects. Our original motivation, however, comes from the studies of 
nonequihbrium steady states. The main difficulty with nonequilibrium steady states is that 
their properties depend not only on the interactions but also on the details of the underlying 
dynamics. In order to get a handle on the effect of dynamics, there has been several works 
|TT|-|n| in which the goal was to determine the effective interactions which are generated by 
the interplay of various dynamical processes. A way to achieve this goal is to investigate 
the phase transitions occurring in nonequilibrium steady states and deduce the effective 
interactions from the universality class the transitions fall into. For example, a combination 
of spin-flips and Levy-flight spin exchanges in a = 1 kinetic Ising model produces a steady 
state with effective long-range interactions and this steady state displays a critical phase 
transition. In this case, the large-distance behavior of the interactions J(r) ~ f^^^" can be 



deduced [jT2[ from measuring critical exponents. Of course, this can be done only if a is in 
a range < cr < 1 where either the thermodynamic exponents or the correlation exponent 
7] depend on a. In order to extend this method to non-integrable interactions (a < 0), we 
must investigate the correlations since measuring thermodynamic exponents yield mean-fleld 
values independently of a. Thus, we can see that in order to extract information about non- 
integrable effective interactions in nonequilibrium steady states, we must flrst understand 
the equilibrium correlations in systems with non-integrable interactions. 

Since mean-fleld theory is expected to describe the critical behavior correctly for a < 0, 
we used it to calculate the correlations for flnite size systems (section Our main result 
is that the singular part of the correlations appears to have universal features which make 
it useful in deducing effective interactions in nonequilibrium steady states. We check our 
mean-fleld results against simulations in section |T| and a flnal discussion is given in section 

m 
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II. CALCULATION OF THE CORRELATION FUNCTION IN MEAN-FIELD 

THEORY 



We consider a one-dimensional Ising model of L = 2N + 1 spins si = ±1 at sites 
i — —N, —N + 1, ...,N and impose periodic boundary conditions. The energy of the system 
is given by 



i<j 



(2) 



and the interaction is assumed to decay with distance as a power law but it is cut off at 

\i-j\^N: 



J/\i-j\^+'' for 1 < |i - il < 
otherwise. 



(3) 



Here J > and the parameter a is assumed to be in the range — 1 < cr < 0. The inhomoge- 
neous external field, Hi is included in order to calculate the susceptibility Xi,j = d{si)/dHj 
which, at temperature = ksT, is related to the two-spin correlation function Cjj = {siSj) 
through 



(4) 



where {Hi} — > means that the susceptibility and, consequently, the correlation functions 
are evaluated in the limit of vanishing fields. 

In the mean-field approximation, Xi,j is obtained from the self- consistency equation for 
the average magnetization: 



(si)=tanh ^(3Ji,i{si)+(3Hi 



(5) 



Taking a derivative of both sides of this equation with respect to Hj and letting {Hi} — > 0, 
we obtain 



(6) 



Equations and from a closed set of equations for Cij and (sj). In zero field, the 

system is expected to be homogeneous and thus Cij = Ci^j and (sj) = m where m is given 

by the stable solution of @. Assuming this homogeneity, equation is solved by Fourier 

transformation. Introducing 

^ 27m 
C{k)= CieM-^kly. k = 7T^-^: n = 0,±l,...,±iV , (7) 



2N + r 



one finds 



^^^^ " 1 - (1 - J)/3J(fc) ' 
where J(k) is the Fourier transform of Jjj- = Ji-j- 

Up to this point, the derivation followed along standard lines. The non-integrability of 
the interaction starts to play role when we try to find the critical temperature, T^. Since 
C{k = is proportional to the fluctuation of magnetization, its divergence determines 
in the thermodynamic limit oo. Thus setting m = in eq.(|D one finds 

^ 27 

ksT, = J(0) = hm 5: ^ . (9) 

The problem with the above expression is that the sum diverges as N^"^ for a < 0. Thus, in 
order to have a thermodynamic limit for the correlation functions, we either have to divide J 
by A^'°"' or have to consider temperatures which are proportional to A^'""'. Since only the ratio 
J/ksT enters the expressions, the two routes are equivalent. Introducing 9 = Tc{N)/T and 
^ = (1 — im?)9^ where Tc{N) is defined through the sum in without taking the N ^ oo 
limit, we obtain the spatial correlation function from the inverse Fourier transform of (^: 
^ l-m^ r 1 JX. cosJ^Z ] ^ ^ 

where Qn^'^) is given by 

=i:^^#=/i:-:^ ■ (11) 

i=i J j=i J 

One can see that correlations in the low- and high-temperature phases are simply related 
by rescaling the temperature and the amplitude by 1 — |jl5|. For this reason we shall 



consider only the high-temperature phase and the critical point and thus set m = and 

e = e. 

Equations (0) and (|1^) can be easily evaluated numerically and the mean-field results 
can be compared to MC simulations (see next Section). First, however, we shall discuss 
the analytic properties of Ci in the "thermodynamic limit" where ^ oo. As can be 



seen from ( [ToD the average of Ci over / is given by [(2A^ + 1)(1 — 9)Y^ and is negligible 
for 9^1. The part of Ci which remains after subtracting the average also diminishes but 
with a smaller exponent. It is of the order of A^"^ and thus no spatial correlations remain 
in the N ^ oo limit for a < 0. It turns out, however, that the correlation function has 
an interesting singular behavior for small x = 1/{2N + 1) which can be well observed in 
finite-size systems. We shall investigate this spatial dependence by amplifying it through 
subtracting the average part of Ci and multiplying the result by 2N + 1: 

The function QN{n) is finite and can be evaluated exactly in the N —>■ oo limit 

7r-r(l-a)cosf _ 
lim QN{n) = n = —n y (13) 



where T{z) denotes the gamma function. Thus we can take the upper limit of the sum in ([12[) 
to infinity and find a scaling function, $(x), which is independent of A^. The singular part of 
$(x) in the limit x — > is calculated by first noting that the identity 2 Y,n=i ^'^^ {27mx) = — 1 
is valid for any x that can be written as 1/{2N + 1) where / is an integer not divisible by 
2A^+1. Using this identity, the sum in (|12D can be rewritten so that it contains the expression 
(1 — 6Da/n^"^)^^ — 1. This expression is then expanded in 6 and, finally, the sums over n 
are changed into integrals. The number of singular terms obtained in this way depends on 
a and, for (j + 1)^^ < |(t| < j^-*^, where j is an integer, we find: 

+ . (14) 

The coefficients A„((t) can be calculated explicitly. For example, the coefficient of the 
strongest singularity, Ai, is given by 



Ai((t) =2-^-"|a| . (15) 

If |cr| = 1/j then the last term in (0) changes from power-law singularity into a term 
of logarithmic singularity, Aj In (x). For example, the case cr = —1/2 yields the following 
singular part 

6 9^ 

$,„,(x) = -= + -ln(x) . (16) 
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One can compare ^n{x) calculated for a = —0.5 and 6 = 0.5 from the exact sum with 



e.g. 2A^ + 1 = 801 and ^singix) obtained from (|T^). One then finds that <l>Ar(x) - $ 



smg 



is constant {A = —2.35) to an accuracy of the order of 1%. Thus the singular part plus 
a constant represents the exact solution to a very good accuracy over the entire range, 
< X < 0.5, which is relevant physically. We checked the "singular-terms-plus-constant" 
description for other values of a and found that the accuracy decreases as a decreases. For 
example, in order to achieve a 1% accuracy at cr = —0.75 and = 1, one must also include 
an additional term and write $(x) ~ Aix~^^^{l + Bx) + C . 

Next, we turn to the calculation of $(x) at the critical point {9 = 1). In systems with 
short-range interactions, the decay of correlations changes from exponential to power law 
as the critical point is reached. Systems with long-range, integrable interactions (a > 0) 
have power-law correlations already in the high temperature phase and the change that 
occurs in these systems when is approached is that the exponent of the power law decay 
changes (becomes smaller at Tc). For example, in the exactly solvable spherical model ||16|| , 



Ci ~ for T > Tc and Ci ~ at Tc. As we shall see below, the non-integrable 

systems are different from the cases discussed above in the sense that the singular structure 
of Ci does not change when 6' = 1 is reached. The scaling function ^sing{x) at criticality is 



obtained by just setting = 1 in eq.(|14|). 



The first question that arises when the 6 = 1 case is considered is the infinity of the 
1/(1 — ^^) term which has been subtracted from (2A^ -|- l)Ci in the definition of $ (see (|12|)). 
Since this term is the sum of the correlation function, it is equal to the magnetization 
fluctuation: 



T3^ = c^(* = o) = l}vTr • (1^' 

where M = J^i^i is the total magnetization of the system. Thus, at the critical point, 
(1 — 6)"^ does not diverge in a finite-size system, it just changes from a quantity of 0(1) 
to 0{N^) where p = 1/2 in the mean- field theory. The apparent singularity of 1/(1 — ^) is 
the result of effectively linearizing eq.(^) when using eq.(^ for calculating the correlation 
function. For our purpose of calculating <l>, the real or apparent singularity of C(A; = 0) 
does not pose a difficulty since $ is defined after subtracting C{k = 0) from the correlation 
function. Note that $(x) is well defined physically: it is the correlation function after the 
average has been subtracted. 

The second question that has to be answered for 6' = 1 is whether any of the terms in 
the sum (|T2|) defining $ becomes singular in the 6^1 limit (After all 1/(1 — was a term 
in this sum before being treated separately and, in general, one may expect long-wavelength 
singularities at the critical point). The answer to this question is no and the reason for 
the absence of singularities is that long-range interactions develop a gap in the "fiuctuation 
spectrum" of the system. In order to see how this happen and why non-integrable interac- 
tions are different from the short-range ones, let us examine the denominator of the terms 
in the sum in eg . ([lOD . For short-range, e.g. nearest-neighbor, interactions one finds in the 
long-wavelength limit 

2-71 / 2TTri \ ^ 

l-eQNin) = l-ecos n^l~e + -i ) . (18) 

' 2iV+ 1 2 V2iV+ 1/ ^ ' 

One can see that the above expression is nonzero for 9 < 1 while it goes to zero at ^ = 1 as 
the square of the wavenumber (for oo). This is the origin of the exponential (power 

law) decay of large- distance correlations at < 1 [9 = 1). 

The structure of the denominator changes for non-integrable potentials (cr < 0). It 
becomes independent of for oo and, using (|13D, one finds 

i-»Q„(„)«i-o^:ni^^)f^ . (19) 
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since D^j = 7r°"r(l — cr) cos ^ < 1 for — 1 < o" < {D^ 1 only for o" — 0), the denominator 
in (p!9D is nonzero even for the smallest n = 1 at 6 < 1. Thus, there is a gap in the spectrum of 
magnetization fluctuations for all wavenumbers, and the singularity of the scaling function 
$(2;) at small arguments (small distances) comes from the non-analytic structure of the 
denominator at large wave-numbers. It should be noted that the above argument also holds 
in the low-temperature phase after replacing ^ by ^ and using the fact that < ^ < 1 ||T5|. 
It is clear from the above discussion that 6 = 1 can be set in the sum ([T^) for evaluating 



^(x), and it is also clear that the result for the singular part ^sing is given by eq.(|Tj) with 
6 = 1. One can again compare $smg and $, the latter evaluated for large but finite systems 
and one again finds 1% accuracy with the ^sing + const description at a = —0.5. The trends 
with the accuracy of the singular-part description is also the same for ^ = 1 as in the case 
of 6* < 1 discussed above. 



III. MONTE CARLO SIMULATIONS 

In this section, we describe the results of addressing the following question: Do the 
results of the mean-field calculation remain valid when fluctuations are included? Since 
no exact solutions are available, we studied the problem by performing MC simulations on 
the long-range Ising model described at the beginning of Section II. We considered cases of 
a = —0.50 and —0.25 at both 6 = 0.5 and 1.0. Single-spin-fiip dynamics with Metropolis 
flip rates was used and we studied systems of sizes 2N + 1 = 201, 401, 801 and 1601. After 
an initial estimate of the relaxation times, we let the system equilibrate and the equilibrium 
correlations were calculated and analyzed. 

Apart from the usual difficulties of simulating a system with long-range interactions 
(calculating the energy changes is time-consuming and no cluster algorithms can be used 
to reduce critical slowing down), there is one extra problem when a cr < system is simu- 
lated. Namely, the definition of the 6 = 1 critical point is not quite obvious. The mean-field 
T^^^{N) defined in (^ has the right asymptotics for — > 00 but for a finite A^, Tc has 
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corrections due to fluctuations. For a sequence of systems of increasing size, one can, in 
principle, use any sequence of Tc(iV) which is approaching T^^{N) in the thermodynamic 
limit. In practice, we defined Tc{N) as the temperature where the macroscopic magneti- 
zation fluctuations have the large- A^" behavior that follows from mean-field theory, namely 
{M^)/{2N + 1) « 1.17..(2A^ + 1)V2. Since {M^)/{2N + 1) = 1/(1 - 6), the above choice 
means that we approach the 9c = 1 critical point as 9c{N) = 1 — [1.17..(2A'" -|- 1)]~^/^. The 
same correction can also be made in the mean-field calculation. This correction, however, is 
small and would not be visible in the Figures discussed below. 

Fig.l displays the raw data for the correlation function at cr = —0.5 and 9c{N). After 
subtracting the average correlation and multiplying the result by 2N + 1 (i.e. calculating 
$) we obtain a collapse of data (Fig.2) if $ is plotted against x — 1/{2N + 1). A blow-up of 
the small x region is shown in Fig. 3. 

In Fig. 4, the MC data is compared with $(x) calculated from the mean-field theory of 
Sec. II, and we observe excellent agreement between simulations and theory. On Fig.5, we 
can see similar agreement away from the critical point {9 = 0.5) for the case of a = —0.5. 
We have also performed simulations with a — —0.25 at both 9 — 0.5 and 9c{N) and we again 
observed both the scaling of the data when ^{x) vs. x was plotted as well as the agreement 
between MC and mean-field results. Thus we conclude that the mean-field theory presented 
in Sec.II correctly describes the correlations of Ising models in which the interactions are of 
long-range, power law form Ji ~ with — 1 < cr < 0. 

IV. FINAL REMARKS 

The main conclusions of our work are that (i) the mean-field theory provides us with 
explicit expressions for the correlation function in case of nonintegrable interactions and (ii) 
the singularity of the correlations for a given a has the same structure at all temperatures. 
The results appear to be not affected by fiuctuations as indicated by Monte Carlo simulations 
for particular values of a and T. 
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We expect that the mean-field theory is also correct for higher dimensional systems where 
the fiuctuations play even less role. This problem is currently being investigated. 
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List of figures: 



Figure 1: 

Raw data from the MC simulation of an Ising chain of length 2N + 1 at the critical 
temperature T = with the long-range parameter set to a = —0.5. The correlation 
function, Ci is plotted against /. 

Figure 2: 

Scaling of data in Fig.l after subtracting the average of Ci over / and multiplying the 
results by 2N + 1, i.e. plotting $Ar(x) defined in eq.(0) against x = 1/{2N + 1). For clarity, 
only 40 points are included from each data sets. 

Figure 3: 

Blow-up of the small x region of Fig.2. All data points in the given x interval are 
included. 

Figure 4: 

Comparison of the scaled MC data (Fig.2) and the mean-field results for the 2N -|- 1 = 
1601 data. 

Figure 5: 

Comparison of the scaled MC data and the mean-field results for a = —0.5 and T = 2Tc. 
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